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MULTIPLE DISCRIMINANTS AND CRITICAL 
VALUES OF A MULTIVARIATE POLYNOMIAL. 


Ruslan Sharipov 


Abstract. A critical value of a function is the value of this function at one of its crit¬ 
ical points. Each critical point of a differentiable multivariate function is described 
by the equations which consist in equating to zero all of its partial derivatives. How¬ 
ever, in general case there is no equation for the corresponding critical value. The 
case of polynomials is different. In the present paper an equation for critical values 
of a polynomial is derived. 


1. Introduction. 

Let f(xi ,... , x n ) be a smooth real valued multivariate function in R” or in some 
open domain C R ra . Critical points of the function / are determined by solving 
the following system of equations with respect to the variables X \, ... , x n : 

df(x i,... ,x n ) _ Q 


dfjxi,... ,x n ) = 

dx n 

They can be maxima, minima, saddle points, and other types of critical points. 

Let v be the value of f(x i,... ,x n ) at one of its critical points given by the 
equations (1.1). Then we can extend the system of equations (1.1) as follows: 


' df(x-i,... ,x n ) _ 

dxi ’ 

< df(x 1 ,... ,x n ) (1-2) 

dx n 

. f(x 1 ,... ,x n ) - v = 0. 


The parameter v in (1.2) is treated as a new variable independent of n ... , x n , 
i. e. (1.2) is a system of n + 1 equations for n + 1 variables. Theoretically, one 
can eliminate the variables x\ ... , x n from the system (1.2) thus producing one 
equation for one variable v. which is the critical value of /: 


F(v) = 0. (1.3) 
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In practice it is rather difficult to derive such an equation. In the present paper 
we derive the equation of the form (1.3) for the case where ,x n ) is a 

multivariate polynomial. For the sake of simplicity, from now on we treat X\ ... , x n 
as complex variables so that they constitute a point in C n . 

2. The case of a univariate polynomial. 

Let’s begin with the case where n = 1. Denoting x\ = x for the sake of simplicity, 
assume that f{x) is a univariate polynomial of n-th degree: 

f(x) = a 0 x n + a\ x 71-1 + ... + a n -i x + a n . 

Then p(x ) = /( x) — v is also a univariate polynomial of n-th degree: 

p(x) = a 0 x n + ai x n ~ x + ... + a n -1 x + (a n - v). (2.1) 

Using p(x), in this case we can write the system of equations (1.2) as follows: 

f p'(x) = 0, 

1 p{x) = 0. 

It is known that a univariate polynomial p{x) vanishes along with its first derivative 
p'(x) at some point x £ C if and only its discriminant is zero (see [1]): 

D p = 0. (2.2) 

The coefficients of the polynomial p(x) in (2.1) depend on v. Therefore the dis¬ 
criminant D p in (2.2) depends on v. And we write (2.2) as 

D p {v) = 0. (2.3) 

The equation (2.3) is a required equation of the form (1.3). 

Theorem 2.1. A complex number v is a critical value of a univariate polynomial 
f(x) if and only if it is a root of the equation (2.3), where D p {v) is the discriminant 
of the polynomial p{x) = f(x) — v. 

Theorem 2.1 follows immediately from the above mentioned basic property of 
the discriminant of a univariate polynomial. 

3. The case of a multivariate polynomial. 

This case is somewhat similar to the previous one. Here we can replace the 
multivariate polynomial /( X\,... ,x n ) with the polynomial 

p(xi, ... , x n ) = f(x i,... ,x„)-v. (3.1) 

Due to (3.1) we can write the equations (1.2) as follows: 

' 9p(x i,... ,x n ) _ 
dx\ 

dp(x i,... ,x n ) = 
dx n 

. p{Xl, ... , Xn) = 0. 




(3.2) 
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The concept of the discriminant of a multivariate polynomial is not commonly 
known. Its definition can be found in [2] (see also [3] and [4]). 

Definition 3.1. The discriminant D p of a multivariate polynomial p(x i,... ,x n ) 
is a polynomial function of the coefficients of p such that its own coefficients are 
integer, such that it is irreducible over Z, and such that D p = 0 if and only if the 
system of equations (3.2) has an least one solution in complex numbers x\, ... , x n . 

As it was said in [2] with the reference to [5], discriminants of multivariate 
polynomials were first considered by G. Boole. The reference [5] is taken from [2] 
“as is”. It looks quite uncertain, since no information on publishers is provided (see 
more detailed historical research in [6] and [7]). 

The coefficients of the polynomial p(x i,... , x n ) in (3.1) depend on v. Therefore 
the discriminant D p depends on v. Hence we can write the equation 

D p (v) = 0. (3.3) 

The equation (3.3) is of the form (1.3). It solves our problem of reducing (1.2) to 
a single equation for v through the following theorem which is immediate from the 
above Definition 3.1. 

Theorem 3.1. A complex number v is a critical value of a multivariate polynomial 
f{x i,... ,x n ) if and only if it is a root of the equation (3.3), where D p (v ) is the 
discriminant of the polynomial (3.1). 

However, the matter is that there is no simple formula for the discriminant D p 
in the case of a multivariate polynomial p. For this reason below we consider a 
different approach to deriving an equation of the form (1.3) from (1.2). 

4. Critical points of discriminants. 

Let p(x, y) be a univariate polynomial of n-th degree in x, i. e. 

p(x, y ) = a n (y) x n + a n _i (y) x n ~ l + ■■■+a 1 (y)x + a 0 (y), (4.1) 

whose coefficients are smooth function of a complex variable y £ C. Assume that 
at some point ( xo,yo) G C 2 , where a n (y o) ^ 0, the following equations are fulfilled: 


p' x ( x o,yo) = 0, p' v (xo,yo) = 0, p(x o ,y o )=0. (4.2) 

Theorem 4.1. If p is a univariate polynomial of the form (4.1) whose coefficients 
are smooth function of a complex variable y and if the equations (4.2) are fulfilled 
at some point (:ro,2/o) G C 2 , where a n (yo) 0, then they imply the equations 


D'(y 0 ) = 0, 


D(y 0 ) = 0 (4.3) 


for the discriminant D of p which are fulfilled at the point yo £ C, where yo is the 
second coordinate of the corresponding point (xo,yo) G C 2 . 

Since a n (yo) 0 in Theorem 4.1, we can divide the polynomial (4.1) by a n (y) 
and proceed to the following monic polynomial: 


P(x,y) = x n + 


a n -i(y) 

a n {y) 


| «i (y) x ao(y) 

a n {y) a n (y)' 


, 71—1 


(4.4) 
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The discriminants of the polynomials (4.1) and (4.4) are related as follows: 

D p (y) = (a n (y)) 2n ~ 2 D p (y). (4.5) 

Relying on (4.5), we can write the monic polynomial 

p(x, y) = x n + a n -i(y) x” _1 + ... + ai(y) x + a 0 (y) (4.6) 

and then reformulate Theorem 4.1 in the following way. 

Theorem 4.2. If p is a monic univariate polynomial of the form (4.6) whose 
coefficients are smooth function of a complex variable y and if the equations (4.2) 
are fulfilled at some point ( Xo,yo ) € C 2 , then they imply the equations (4.3) for the 
discriminant D of p which are fulfilled at the point yo £ C, where yo is the second 
coordinate of the corresponding point (xo,yo) £ C 2 . 

Theorems 4.1 and 4.2 are equivalent to each other due to the relationship (4.5). 

5. Proof of Theorem 4.2 in the case of a quadratic polynomial. 

Assume that n = 2 in (4.6). Then p(x,y) is a quadratic polynomial: 

p(x,y) = x 2 +ai(y)x + a 0 {y). (5.1) 

One can easily calculate the partial derivatives of p(x, y): 


P'x(x,y) = 2x + a 1 {y), p' y (x,y) = a[{y) x + a' Q {y). (5.2) 

Substituting (5.1) and (5.2) into (4.2), we derive a system of three equations: 


2x 0 + ai(yo) = 0 , 

al(yo)x 0 + ao(y 0 ) = 0, 

(x 0 ) 2 + ai(yo) xq +a 0 (y o) = 0. 

The discriminant of the polynomial (5.1) is calculated as follows: 

D{y) = (ai(y)) 2 -4a 0 (y). 

Its derivative D'(y) is also easily calculated: 

D'(y) = 2oi (y) a[(y) - Aa' 0 {y). 

Note that the first equation in (5.3) can be resolved with respect to xo- 

ai(yo) 

*o =-—• 

Substituting (5.6) into the third equation (5.3) we derive 

_(fLMi + oo(w) = „. 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 
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Comparing (5.7) with (5.4), we see that (5.3) implies 


D(yo) = 0. 


(5.8) 


Then we substitute (5.6) into the second equation (5.3). As a result we get 


aijyo) a[(yo) 
2 


+ a-oiyo) = 0 . 


Comparing (5.9) with (5.5), we see that (5.3) implies 


D'(y 0 ) = 0. 


(5.9) 


(5.10) 


The rest is to note that (5.10) and (5.8) coincide with (4.3) and conclude that (4.2) 
implies (4.3) in the case of the polynomial (5.1). Thus, Theorem 4.2 is proved for 
the case of quadratic polynomials. 


6. Proof of Theorem 4.2 in the case of a double root polynomial. 

If y = yo is fixed, the equations (4.2) mean that the polynomial p vanishes along 
with its derivative p' x at the point x = Xo, i- e. it has a multiple root at this point. 
The case of a double root is the most simple in this situation. 

Let x\ ... , x n be roots of the polynomial (4.6). They depend on y, i. e. 

Xi = Xi(y ), i = 1, ... , n. (6.1) 

It is known that roots of a univariate polynomial are continuous functions of its 
coefficients (see [8]). Hence in our case the roots (6.1) are continuous functions of 
y. Two of them tend to Xq as y —> yo- Without loss of generality we can set 

xi(y)x 0 and x 2 (y) ^ x 0 as y -A y 0 . (6.2) 

Assume that the roots £ 3 (yo), ■ ■ ■ , x n (yo ) are distinct and different from the double 
root xi(yo) = x 2 (yo) = xq. Under this assumption the roots 


X3(y), ■ ■ ■ , x n (y) (6.3) 

are smooth functions of y in some neighborhood of the point y = yo (see [8]). Unlike 
them the roots X\ (y) and x 2 (y) in (6.2) are not necessarily smooth, though they 
are continuous functions of y in this neighborhood of the point y = yo- 

Using the roots (6.2) and (6.3), we define two polynomials which are two com¬ 
plementary factors of the initial polynomial p(x,y) in (4.6): 

2 n 

<l(.x,y) = Y\(x-Xi{y)), r(x,y) = JJ(x - x*(y)). (6.4) 

i =1 i=3 

Indeed, from (6.4) we easily derive the equality 


P(x,y) = q{x, y) r(x, y). 


(6.5) 
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The coefficients of the polynomial r(x , y ) in (6.4) are smooth functions of y since 
they are elementary symmetric functions of the smooth roots (6.3) (see [1]). The 
polynomial q(x, y) in (6.4) is quadratic: 

q{x , y) = x 2 + /3i (y) x + f3 0 (y ). (6.6) 

Its coefficients are also smooth functions of y since due to (6.5) the polynomial (6.6) 
can be produced by means of the polynomial long division algorithm (see [9]) with 
the monic polynomial r(x, y) as a divisor: 

q(x,y) = p(x,y) + r{x,y). 

Let’s recall that the roots £3(2/0), ■ • ■ , x n(yo) are different from the double root 
xq. Therefore from (6.4) we derive the inequality 

n 

r(x 0 ,yo) = JJOeo - Xi(y 0 )) ^ 0. 

2—3 

Similarly, applying (6.2) to q(x,y) in (6.4), we derive the equality 

q{xo,yo) = 0 . 

Differentiating the equality (6.5), we find that 

p' x ( x o,yo) = qx(x 0 ,y 0 )r{x 0 ,yo) + q(x 0 ,y 0 ) r' x {x 0 ,y 0 ), 

Py{xo,yo) = qy{xo,yo)r(x 0 ,y 0 ) + q(xo,yo)ry{x 0l yo)- 

Now, applying (6.8) to (6.9), we obtain the equalities 

p' x ( x o,yo ) = qL( x o,yo)r{x 0 ,yo), 

Py{x o, 2/o) = q' y {x o, 2 /o) r(x Q , y 0 ). 

If the equations (4.2) are fulfilled, then, using (6.7), from (6.10) we derive 

q'x{ x o, 2/o) = 0, gd(£ 0 ,2/o) = o. (6.11) 

The equalities (6.11) combined with (6.8) mean that the equations (4.2) for p, 
once they are fulfilled, imply similar equations 

qL( x o,Vo) = 0, qy{x o ,y o )=0, q(x o ,y o ) = 0 (6.12) 

for the quadratic polynomial (6.6) whose coefficients are smooth functions of y. 
Theorem 4.2 is already proved for quadratic polynomials. Therefore, applying this 
theorem to q{x,y) and (6.12), we derive the equalities 


(6.7) 


( 6 . 8 ) 


(6.9) 


D'q{yo) = 0, 


D q (y 0 ) = 0, (6.13) 


where D q (y) is the discriminant of the polynomial q(x,y) in (6.6). 
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Let’s proceed to the discriminant D(y) of the polynomial p(x,y) in (4.6). In 
terms of its roots the discriminant D(y) is given by the formula 


D (y) = JlMv) - x Av)f ( 6 - 14 ) 

i<3 

(see [1]). If we denote through D r (y) the discriminant of the polynomial r(x,y ) in 
(6.4), then we can factorize the discriminant (6.14) as follows: 

2 n 

D (y) = D q(y) D r(y) n ~ x Ay)f- ( 6 - 15 ) 

i =1 j =3 


Taking into account the formula for q(x,y) in (6.4), we can transform (6.15) as 


D{y) = D q(y) D r(y){Y[q(xj(y),y)j . 

kj=3 ' 

Relying on (6.16), lets introduce the following notation: 


(6.16) 


a{y) =D r {y) ( Y[q(xj(y),y) j . 

^ j=3 ' 


(6.17) 


The function a(y) in (6.17) is a smooth function of y since the coefficients of the 
polynomial q(x,y ) in ( 6 . 6 ) and the roots *3 (y), ... , x n (y) of the polynomial r(x,y ) 
in (6.3) are smooth functions of y. Applying (6.17) to (6.16), we get 


D(y) = D q (y)a(y). (6.18) 

Both multiplicands D q (y) and a(y) in (6.18) are smooth functions of y. Therefore, 
differentiating (6.18), we derive the following formula: 

D \yo) = D q(yo ) a(y 0 ) + D q (y 0 ) a'(y 0 ). (6.19) 

The rest is to apply (6.13) to (6.18) and (6.19) and derive (4.3). Thus we have 
proved that (4.2) implies (4.3) in our present case. Theorem 4.2 is proved in the 
case of a polynomial (4.6) with exactly one double root and the other simple roots. 


7. Proof of Theorem 4.2 in general case. 

Note that the polynomial p(x, y) in (4.6) depends on y through its coefficients. 
Therefore the derivative D'{y) of its discriminant is calculated is follows: 


n } on 

C '(») = £gV “.'W 


2=0 


(7.1) 


Let’s denote Sp(x,y) =p' y {x,y). Then, differentiating (4.6), we get 


Sp(x, y) = a,n_i(y) x n 1 + ... + a[(y) x + a' 0 (y). 


(7.2) 
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Let’s denote a; = a[{yo), h = a[{yo) and SD = D'(yo). Then, substituting y = yo 
into (4.6), (7.2) and (7.1), we obtain two polynomials 

p(x) = x n + a n -1 £ n-1 + ... + <zi x + a 0 , 

Sp(x) = b n - 1 x n ~ x + ... + bi x + b 0 


with purely numeric coefficients and the purely numeric quantity 

dD 

dai 




2=0 


(7.4) 


associated with the polynomials (7.3). The following definition is terminological. 
It is designed for the sake of beauty. 

Definition 7.1. The polynomial 5p(x) in (7.3) is called a first variation of the 
polynomial p{x), while the numeric quantity (7.4) is called the first variation of the 
discriminant D of p associated with 5p(x). 

It turns out that the functional nature of the coefficients of p(x , y) in (4.6) is 
inessential in Theorem 4.2. This theorem can be reformulated as follows. 

Theorem 7.1. If a monic polynomial p(x) in (7.3) has a multiple root of mul¬ 
tiplicity m ^ 2 and if it shares this root with its first variation 8p(x), then the 
discriminant D of p(x) vanishes along with its first variation 5D in (7.4). 

Theorem 7.1 is equivalent to Theorem 4.2. The equivalence can be established 
using the linear functions cq(y) = ai + bi (y — yo) in (4.6). The result of the previous 
section means that we have already proved Theorem 7.1 for any monic polynomial 
p(x) with exactly one double root and the other simple roots. 

Let p{x) be a monic polynomial with at least one multiple root x$. Then it can 
be produced as a limit of some sequence of monic polynomials p s {x) with exactly 
one double root x = xo and the other simple roots. Assume that p(x) shares its 
multiple root x = Xq with its first variation 8p(x) in (7.3). Then 


Sp(x 0 ) = 0. 


(7.5) 


The equality (7.5) is written as a linear relationships with respect to the coefficients 
bo, , b n - 1 of the polynomial 8p(x) in (7.3): 

Xq 1 b n -1 + ... + Xo &i + b 0 = 0. (7-6) 

The coefficients of the linear combination (7.6) depend only on the root xo, which 
is common for p(x) and for any polynomial in the sequence p s (x). This means that 
p{x) and the polynomials p s (x) share the same first variation 5p(x) obeying the 
relationship (7.5). This first variation Sp(x) combined with each p s (x) according 
to (7.4) produces a numeric sequence 5D S . The discriminant D and its partial 
derivatives in (7.4) are smooth functions of ao, ■ ■ ■ , On-i- Therefore we have 

SD = lim SD s . (7.7) 

s—too 

The rest is to apply Theorem 7.1 to each p s {x) combined with Sp(x). This yields 
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SD s = 0. Substituting 5D S = 0 into (7.7), we derive the required result 5D = 0. 
Thus, Theorem 7.1 is proved in its full generality. The same is true for Theorems 4.2 
and 4.1, which are equivalent to Theorem 7.1. 

8. Application to critical values. 

Having been equipped with Theorem 4.1, now we return to our initial problem 
of deriving an equation of the form (1.3) from the equations (1.2) in the case 
of a multivariate polynomial f(x i,... ,x n )- Passing from f{x i,... ,x n ) to the 
polynomialp(xi,... , x n ) in (3.1), we treat p(x i,... , x n ) as a univariate polynomial 
with respect to the variable x = x n and treat the other variables X\. ... , x n -\ and 
v as parameters. Then we calculate the univariate discriminant D of the polynomial 
p{x i,... , x n ) with respect to the last variable x = x n \ 

P(x i,... ,x n -i) = D(x n ,p(x i,... ,x n )). (8.1) 

The univariate discriminant D in (8.1) acts as a nonlinear operator sending the n- 
variate polynomial p(x i,... , x n ) to the (n — l)-variate polynomial p(x i,... , x n _i). 
Setting y = Xi and applying Theorem 4.1 for each i = 1, ...,n — 1, we derive the 
following system of equations from the equations (3.2): 

' dp{x i,... ,x n -i) _ 
dx\ 

< dp(x i,... , x n -\) _ (8.2) 

dXn-X 

, p{x 1 ,... , x n -\) = 0. 

The structure of the equations (8.2) is the same as the structure of the equations 
(3.2). Therefore we can apply the operator (8.1) repeatedly: 

DD P = D(x!,D(x 2 ,. ■ • ,D(x n ,p(x i,... ,x n ))...)). (8.3) 

Definition 8.1. The numeric quantity DD p introduced through the formula (8.3) 
is called the multiple discriminant of a multivariate polynomial p(x ±,... , x n ). 

Generally speaking, the multiple discriminant DD p is different from the discrim¬ 
inant D p of a multivariate polynomial introduced in Definition 3.1. Unlike D p , the 
formula (8.3) provides a clear algorithm for calculating DD p . 

The coefficients of the polynomial p(x i,... , x n ) in (3.1) depend on v. Therefore 
the multiple discriminant DD p depends on v. Hence we can write the equation 

DD p (v) = 0, (8.4) 

where DD p (v) = D(xi, D(x 2 , ■ ■. ,D(x n ,f(x 1 ,... ,x n ) -v))...)). 

Theorem 8.1. If a complex number v is a critical value of a multivariate polyno¬ 
mial f{x i,... ,x n ) : then it is a root of the equation (8.4). 

Theorem 8.1 is proved by applying Theorem 4.1 repeatedly. Theorem 8.1 is 
similar to Theorem 3.1, but it is somewhat weaker, i. e. each critical value of a 
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multivariate polynomial f(x i,... ,x n ) is a root of the equation (8.4), but not each 
root of the equation (8.4) is a critical value of f(x i,... , x n ). 

9. Conclusions. 

Theorem 8.1 along with the equation (8.4) and the formula (8.3) constitutes 
the main result of the present paper. It can be applied in testing positivity of 
multivariate quartic forms and forms of higher degrees. 
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